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pure two-dimensional quantum gravity with an interaction. It permits us to study the 
^ I intermediate regimes between the gravitating and flat metric. Flat space is modeled by 
. 5^ I a regular square lattice, while localised curvature is introduced through lattice defects. 

No "flattening" phase transition is found with respect to the coupling: the infrared 
behaviour of the system is that of pure gravity for any finite R^ coupling. In the limit 
of infinite coupling, we are able to extract a scaling function interpolating between pure 
gravity and a dilute gas of curvature defects on a fiat background. We introduce and 
explain some novel techniques concerning our method of large N character expansions and 
the calculation of Schur characters on big Young tableaux. 
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1. Introduction and overview 

Two-dimensional random geometry is now placed at the heart of many models of 
modern physics, from string theory and two-dimensional quantum gravity, attempting 
to describe fundamental interactions, to membranes and interface fluctuations in various 
problems of condensed matter physics. Therefore it is quite important to study in detail the 
basic universal properties of random geometries, such as their fractal nature, their phase 
diagram, critical phenomena and correlations of physical quantities of geometric origin. In 
addition, one might hope to discover entirely new mechanisms for phase transitions. 

Over the last fifteen years, considerable progress has been achieved in the under- 
standing of noncritical strings, or 2D quantum gravity in the presence of matter fields 
with central charge c < 1. It was based on two different approaches, completing and 
justifying each other. 

The first one, based on the continuous treatment of 2D metric fluctuations by means 
of quantum Liouville theory, was introduced by Polyakov 0. 

The second one, based on a discretisation of the two-dimensional metric in terms of 
random planar graphs, was first proposed in [Q, and [If. 

With the help of powerful matrix model techniques this approach allowed for the first 
time the exact calculation of critical exponents in quantum gravity; both for pure 2D 
gravity without matter 0,0, and, subsequently, for various forms of conformal matter 
with c < 1 (e.g. c = —2 and c = ^ §). These results were successfully confirmed and 
supplemented by the continuous approach in 0, and then in [jTOf . 

Following this breakthrough, numerous further developments have led to a rather full 
understanding of the subject. 

Still, we think that some important questions were left unclarified. In this paper, we 
are addressing one of them: what will happen in the case of pure 2D gravity if we introduce 
coupling constants favoring the fiat configurations among the ensemble of 2D metrics? Can 
we achieve a transition to a new phase of essentially fiat metrics? Or shall we always find 
- on sufficiently big distances - the familiar behaviour of 2D gravity, whatever these new 
couplings are? If the last scenario is true, do we have some universality for intermediate 
scales, when the size of the 2D universe is of the order of the flattening scale? 

On first sight, this question appears to have a simple answer, at least from the point 
of view of the continuous theory. The Euclidean path integral is 



Apart from the cosmological term (controlled by /u) and the (topological) Einstein term 
proportional to the curvature Rg, it does not seem to be meaningful to put further terms 




(1.1) 
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into the action of 2D gravity on dimensional grounds. The simplest term one might want to 
consider is -^R^- The bare coupling constant /3o is however dimensionful and thus should 
be proportional to the cutoff squared. So it is small and in principle should be dropped, 
as well as any further higher derivative terms. 

On the other hand, no argument has ever been given to exclude the following alter- 
native scenario: increasing the coupling (/3o — in our notation), the theory might 
entirely restructure nonperturbatively. In this case the Liouville approach, which is known 
to correctly quantise the theory ( p. . 1|) in the absence (i.e. /3o = oo) of higher derivative 
terms, might no longer yield a good description of the theory. Some calculational at- 
tempts to address this question within the framework of Liouville theory are inconclusive 
precisely because they merely perturb by the R^ term in question (see e.g. [^,[0). 

In this paper, we address and resolve some of these questions in a completely non- 
perturbative way. We take as a starting point the by now very familiar discrete definition 
of the path integral (|1.1|), known to correctly quantise the theory in the absence 
of higher derivative terms. We then introduce couplings whose naive continuum limit pre- 
cisely induces the higher terms indicated in ( |1 . 1| ) . Technically, this is done by solving an 
unusual matrix model that generates the dually weighted graphs proposed and studied in 
our earlier papers |^ and 0]. 

The paper is organized in the following way: 

In section 2, the discretised model of dually weighted graphs is formulated and the 
main physical results, including the absence of a "flattening" phase transition, the joint 
scaling function of the R^ coupling and the cosmological constant and their infrared and 
ultraviolet asymptotics are presented and discussed in a non-technical fashion. We physi- 
cally interpret our model as a statistical mechanics system of a "gas" of point-like curvature 
defects surrounded by flat two-dimensional space. Curvature is quantised, i.e. localised on 
defects whose deficit angle equals an integer times some fixed value. 

In section 3, a review of the character expansion approach to the models of dually 
weighted graphs (worked out in and 0) is given. 

In section 4, the exact solution of the present model is derived in some detail. The 
solution is rather complex but nevertheless explicit. Elliptic functions, already crucial in 
0, play a central role. 

Section 5 demonstrates how to extract the physical results from the involved expres- 
sions of section 4. 

Section 6 is devoted to conclusions and open questions. We also discuss the important 
open issue of the universality of our results and point out that fine-tuning the curvature 
weights - an in principle feasible extension of the present work - might result in new phases 
of 2D gravity. 



2 



Finally, appendix 1. contains the calculation of the Schur characters of GL{N) in 
the limit of a big Young tableau while some more technical details concerning the exact 
solution and its verification are given in appendix 2. 

2. Definition of the model and physical results 

In the particular model used here the two-dimensional random geometry is described 
by planar graphs (with the topology of the sphere) built from plaquettes ( "quadrangula- 
tions"). We are summing over all such graphs with the weights t2,t4,tQ, corresponding 
to vertices with 2,4,6,... neighboring plaquettes, respectively. The partition function is 

^=E n t^*"'-^ (2-1) 

G V2g€G 

where V2q are the vertices with coordination number 2q and #f 2g are the numbers of such 
vertices in the given graph G. Note that one cannot sum up these graphs by the usual 
one-matrix model, since we need to control both vertex and face coordination numbers. 
The discrete, curved manifold is thus described by a graph made from flat squares. The 
curvature appears exclusively in the vertices and corresponds to the deficit angles around 
the vertices: 7r{2 — q). Thus two neighboring plaquettes contribute positive curvature, four 
correspond to zero curvature and six and more result in localised negative curvature (see 
Fig 2.1). 



Fig. 2.1 Flat space and curvature defects. 

For technical reasons, it proved useful to choose a particular parametrisation of the 
above weights t2q'- 

t2 = V\t,U = A, = Ai ^, ... t2q = (^)('?-2). (2.2) 
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With these weights, it is easy to prove, using Euler's theorem, that the partition sum ( |2.1j ) 
becomes 

Z{t, A, /3) = y /32(#-2-4)^ (2.3) 



G 



where A is the number of plaquettes of the graph G and the number of positive 
curvature defects. Note that the latter are balanced by a gas of negative curvature defects. 



whose individual probablities are given in (|2.2 



We expect this model to describe pure gravity in a sufficiently large interval of /3, 
after tuning the bare cosmological constant A (controlling the number of plaquettes) to 
some critical value XdP)- On the other hand, for A fixed and /3 = we entirely suppress 
curvature defects except for the four positive defects needed to close the regular lattice 
into a sphere. It is thus clear that P is the precise lattice analog of the bare curvature 
coupling (3o in (|1.1D . The phase (3 = of "almost flat" lattices - very different from pure 
gravity - was discussed in detail in . 

Let us now summarize the main physical results following from the exact solution (for 
general A and (3) of this model: 

1. A long debated question was whether models of the present type undergo a "flat- 
tening" phase transition at a flnite, non-zero critical value of /3 = /3c. The weak coupling 
region (3 > f3c would then correspond to the standard phase of pure gravity while a puta- 
tive novel "smooth" phase of gravity might exist either at f] = (3c or in the entire interval 
< (3 < Pc- This would constitute an existence proof of continuum 2D gravity. We 
flnd, to the contrary, that there is no "flattening" phase transition at non-zero 13. For any 
given (3 we flnd the powerlike scaling of standard pure gravity on large scales. This means 
that no matter how flat the system is on small scales (of the order of f3~ ^ ) , it destabilizes 
in the infrared into the familiar ensemble of highly fractal "baby-universes" . 

2. The dependence of the partition sum ( |2.3| ) on [3 and the lattice cosmological 
constant A in the vicinity of the flat phase /3 ~ and close to A ~ Ac is given by a simple, 
(presumably) universal scaling function /(x) (deflned through Z{t,X,(3) = fi^)) 
reflecting the transition from flat space to pure gravity: 

/(x)=x^-^x4 + ^x2- A-x (x^-l)^ (2.4) 

where the scaling variable x is given, to leading order, by 

x^^i^^. (2.5) 
We can distinguish the following features: 



4 



(a) There is a degree | singularity at x = 1, correctly reproducing the universal 
string susceptibility exponent 7s = — ^ of pure gravity 0,0. In view of eq.( p.5| ), the 
critical value of the lattice cosmological constant A is therefore given to leading order by 
Ac = 1 — + Therefore, in view of (|2.3| ), the characteristic growth of the random 
surfaces as a function of the lattice area A (= number of plaquettes) is given by 

Z{t,A,(3) ~ ile^^^ A-i. (2.6) 

132 

For any non-zero /3 we do have exponential growth of the number of surfaces, but one has 
to go to larger and larger scales (i.e. use more and more plaquettes) to be able to take the 
continuum limit. If f3 is exactly zero there is no longer any exponential growth and no pure 

9 

gravity continuum limit is possible. The prefactor /3~2 is found in the exact calculation 
in section 5; we are not sure whether it is universal. 

(b) We further see that taking (3—^0 before the limit A — Ac corresponds to the 
limit X — ^ oo. In this limit one finds f{x) ~ + 0(-^), that is, the characteristic 
critical behavior of 2D gravity "silently" disappears and we recover a power series in ^ 
corresponding, in view of (p75|), to a perturbative expansion in lattice defects (3. In this 
limit the characteristic growth of surfaces as a function of area A is 

Zit,A,(3) ~ { A + 0{(3^A^) ). (2.7) 

The leading order corresponds precisely to the almost flat lattices (with exactly four pos- 
itive defects) of 0. The corrections are interpreted as insertions of negative defects, 
balanced by further positive defects. The typical shape of the surfaces in this limit is 
a generalisation of the one we found for "almost flat" graphs in [0: Long, fllamentary 
cylinders growing out from every negative curvature defect. 

(c) It is easy to prove that the scaling function (|2.4[ ) is the simplest possible function 
with the limiting properties discussed in (a) and (b). 

The above results might be interpreted in terms of a continuum model of quantised 
curvature defects, in which the localised defects move around like particles in a gas on a flat 
background space. The deficit angle, A9, of a defect can take on the values = tt, and 
—TV. A positive curvature defect is surrounded by a conical geometry, whereas a negative 
curvature defect corresponds to a saddle-type insertion (see Fig. 2.1). The higher order 
negative curvature defects {—2tv and higher) would not be expected to play a role in this 
limit (the entropy from moving two low order defects around would completely dominate 
that from a single higher order defect). The coupling /3 can be interpreted as a fugacity 
controlling the number of defects. The flat space limit (3^0 consists of four defects of 
degree tt moving around with respect to one another. Varying the fugacity, /3, allows one 
to smoothly interpolate between flat space, ( |2.7| ) (with four defects), and pure gravity ( |^ ) 
(with an inflnite number of defects). 
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3. Review of the technique 

The appropriate techniques for deahng with models of duaUy weighted planar graphs 
were pioneered in two previous publications [|l|] , . Let us recall the key steps of the method 
and collect the crucial formulas allowing the exact solution of the model. The method is in 
fact applicable to a situation more general than the one written in eq. ( |2.1D , to be analysed 
in the present work. The partition function of general dually weighted planar graphs G is 
defined to be 

z{t*,t) = J2 n (3-1) 

G v*^,vqeG 

where v* , Vq are the vertices with q neighbours on the original and dual graph, respectively, 
and ^v* , are the numbers of such vertices in the given graph G. This expansion is 
generated by the following matrix model: 

Z{t*,t) = J VM e-^^' + ^ (3.2) 

with 

Vb{MA) = V r TtB'' (MA)''. (3.3) 

k=i 

The matrices A and B are fixed, external matrices encoding the coupling constants through 



= i i- Tr S*? and t„ = 1 i- Tr A^. (3.4) 

q N ^ q N ^ ' 



The model generalises, for A ^ 1, the standard one matrix model first solved by Brezin, 
Itzykson, Parisi and Zuber [^|. It may no longer be successfully treated by changing to 



eigenvalue variables; a reduction to variables is nevertheless possible: One expands the 
potential in terms of the characters on the group: 



N 



v-i cx: 



.H^^TTb;^ - X«(M^). (3.5) 



Here Bi and {MA)j are the eigenvalues of the matrices B and MA. The first step involves 
rewriting the sum over /c as a double sum over all the eigenvalues of the matrices B and MA 
of — ln(l — Bi{MA)j). Exponentiating the log then gives the product in the numerator. 
The second step uses a group theoretic result to write the product in terms of a sum over 
characters. The character is defined by the Weyl formula: 

det. ,Ao!l^) 

X{hM) = ^ (3-6) 
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where the set of {h} are a set of ordered, increasing, non-negative integers, A(a) is the 
Vandermonde determinant, and the sum over R is the sum over aU such sets. The i?'s 
label representations of the group U{N) and the sets of integers {h} are the usual Young 
tableau weights defined hy hi = i — 1 + ^boxes in row i (the index i labels the rows in 
the Young tableau, i = 1 corresponding to the first row) . Note that the restriction on the 
allowed Young tableaux that any row must have at least as many boxes as the row below 
implies that the {hi} are a set of increasing integers: 

hi+i > hi. (3.7) 

Substituting equation ( |3.5|) into the integral in equation ( p.2|) , we can now do the angular 
integration using the key identity 

J {vn)H XRi^Mn^A) = d-^^ xr{m) xr{A), (3.8) 

where da is the dimension of the representation given by da = A(/i)/ fl^j^^ z!), and 
arrive, after performing a Gaussian integral over the eigenvalue degrees of freedom, at the 
expression 

^ r\ (he _ ivi/7°ii 

Zit.n = c Y: n )h^-ho; XwiA) X{Hy{B). (3.9) 

The sum is taken over a subclass of so-called even representations. These are defined as 
possessing an equal number of even weights /if and odd weights h° (since the mentioned 
Gaussian integration vanishes if the latter condition is not satisfied). The formula |15] was 



originally discovered by Itzykson and Di Francesco [jT5| by summing up "fatgraphs" , using 
purely combinatoric and group theoretic arguments. 

In our first paper we demonstrated how to take the large N limit of this expansion. 
In this limit the weights jjhi condense to give a smooth, stationary distribution dh p{h), 
where p{h) is a probability density normalized to one. For simplicity we restrict our 
attention to models in which the matrices A and B are such that traces of all odd powers 
of A and B are zero. This means that the our random surfaces are made from vertices 
and faces with even coordination number. As was discussed in this ensures that the 
support of the density p{h) lies entirely on the real axis and thus simplifies the solution of 
the problem. The matrix A will satisfy this condition if we introduce a ^ x ^ matrix ^/a 
in terms of which A and the character X{h}iA) are given by 



A = 





-x/a 



and X{h}{A) =X{!^}{a)X{h^}{a) sgn[Y[{hi -h°)]. (3.10) 



^,3 
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In order to study the transition from flat to random graphs we choose the potential to be 



Va{MA^) = J2^ TrfA^'^] {MA,f\ (3.11) 
k=i 

Here A4 is deflned to satisfy Tr[(A4)'^] = Ndk,4 and A is as deflned in ( |3.10| ). We are 
thus studying surfaces made up entirely from squares. A weight t2k = Tr[A'^''] is assigned 
whenever 2k squares meet at a vertex (see Fig. 2.1). We are therefore precisely considering 
the situation deflned in eq. (|2.1|) . In order to investigate (|3.9|) in the large limit, one 
attempts to locate the stationary point. This leads to the saddlepoint equation 

2F(h) + r dh' = -\nh. (3.12) 

Jq h-h' 

As has been discussed in detail in |1[], this equation actually does not hold on the entire 
interval [0, a], but only on an interval [6, a] with < 6 < 1 < a: Assuming the equation to 
hold on [0, a] would violate the implicit constraint p{h) < 1 following from the restriction 
( p.7|) . The density is in fact exactly saturated at its maximum value p{h) = 1 on the 
interval [0, a]. The solution of ( |3.12|) requires the knowledge of the large A^ limit of the 
variation of the characters in eq.( |3.1Cl| ): 

F{hk) = 2-^ In \' \ . (3.13) 

A rather general method for its determination has been one of the main technical achieve- 
ments of our previous work. In fact, we found the following simple result: Introduce the 
character function F{h) and the resolvent H{h) as, respectively, the large A^ limit of ( |3.13| ) 
and 

We found the weight resolvent H{h) to be very closely related to the standard matrix 
model eigenvalue resolvent. It provides a direct link between the statistical distribution of 
Young weights and the correlators of the model: 

N q J 2m ^ ^ 

Here the contour encircles the cut of H{h). Further introduce the function G{h) as 

G'(/i) = e^('^)+^('^). (3.16) 
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Its importance stems from the fact that it relates in a simple way the introduced functions 
F{h),H{h) and the coupling constants t2q 



Using this equation as well as an alternative representation of the Weyl character (|3.6| ) 



X{h}{A) = del{Ph,+i-i{e)) (3.18) 

(fc.i) 

as a determinant of Schur polynomials Pn{d) defined through 

oo ^ 

= J2^''Pn{0) with e, = - Tr[A'], (3.19) 

n=0 ^ 

one derives ^ from the "equation of motion" 



d N 

-Pn{e) = Pn-qie) with e, = —t2q (3.20) 



the equation 



^-i-E|f+f:|a^K^m(«))^^ (3.21) 



. , - — ''2nr 

9=1 9=1 ^ 



where the coefficients of the positive powers of G in (|3.21|) are directly related to the 
correlators of the matrix model dual to ( |3.11j ), i.e. the model with potential VA^^iMA) — 
I (MA)': 

Here the negative powers of G are fixed by eq. (|3.17| ) while the positive powers are deter- 
mined through eq.( p.20|) . If we write formula (|3.21|) in the form h Xh{t) = {J2q>o ^ t2q + 
Sq>o 7^ Xh{t) we see that the operator h acts on the character like an operator 
of the derivative of the bosonic field ^ (piG) since the commutation relations of (j){G) 
for different G's are those of a bosonic field0. This fact suggests that the models of du- 
ally weighted planar graphs might have a description in terms of integrable hierarchies of 
differential equations, just as for the (much simpler) case of ordinary matrix models. 

We have also assumed for the moment that only a finite number Q of couplings are 
non-zero (i.e. t2q = for q > Q)- Furthermore, we were able to show in our second work 
P] that (|3.21D implies the functional equation 



e"-' = — n (3.23) 

9=1 



^ We acknowledge several conversations with I. Rostov concerning this point. 
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where the Gq{h) are the first Q branches of the multivalued function G{h) defined through 
( p.21|) which map the point /i = oo to G = 0. The saddlepoint equation (|3.12| ), together 
with (|3.23| ), defines a well-posed Riemann-Hilbert problem. It was solved exactly and in 
explicit detail in for the case Q = 2, where the Riemann-Hilbert problem is succinctly 
written in the form 

2m) + //W^-l.(-ii) ^^^^^ 

2F{h)^${h) = -In/i, 

where ^{h) denotes the real part of H{h) on the cut [6, a] and ^{h) denotes the real part 
of F{h) on a cut [— oo, c] with c < b. This case corresponds, in view of the potential ( |3.11 



to an ensemble of squares being able to meet in groups of four (i.e. fiat points) or two 
(i.e. positive curvature points). We termed the resulting surfaces "almost flat". It turned 
out that all the introduced functions could be found explicitly in terms of elliptic functions. 
We also wrote down the Riemann-Hilbert problem for the more complicated case Q = 3. 
Allowing in addition a non-zero negative curvature coupling te, it already captures the 
transition from flat to random surfaces. The explicit equations are 

2F{h)+:^{h) =-ln/i 

h (3-25) 

Fi{h) + F2{h) + Fs{h) + 2H (h) =-ln{- ' 
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Its exact solution can, in principle, be obtained as well. Some of the steps of the solution of 
( p.25|) as well as for the general case of an arbitrary number of couplings will be presented 
in appendix 1. Unfortunately it involves already functions more general than elliptic, 
making the solution much less explicit and transparent. Luckily, however, we are not 
forced to solve the system ( |3.25D in order to analyse the problem of the transition from 
flat to random lattices. Indeed, we merely need to perturb our almost flat lattices by any 
operator containing negative curvature. This physical observation will be used in the next 
section to achieve the long-sought solution to our problem. 



4. Exact solution of discrete two-dimensional gravity 

In the review of the method of the previous section we recalled that the analytic 
structure of the multivalued function G{h) (|3.16| ) is in part determined by the weights t2q. 
Specifically, the structure of the physical sheet and of all the sheets attached to it by the 
cut of e^^'^^ is the same as that given by (i) dropping the positive powers of G in ( p.21| ) and 
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(ii) inverting this truncated equation to obtain G as a function of h. In this observation 
permitted us to solve the flat space limit. There the equation ( |3.21| ) reads: 



h — 1 = + + positive powers of G. 
G G 



(4.1) 



Dropping the positive powers of G leaves us with a quadratic equation for G. We thus 
deduced that the physical sheet and all the sheets attached to it by the cut of e^^'^^ have a 
quadratic structure. In other words, there is only one sheet attached by the cut of e^*^^-*. 
This simple structure allowed us to solve the flat space limit directly. 

The key observation that allows us to solve the full problem is that we can choose a 
more general set of weights, including those needed to introduce negative curvature defects, 
whilst still preserving this simple two sheet structure. Labeling the flrst two weights as 
in : 7 we now include all the even with q > 3, assigning them the following 
weights: t2q — ^46^"^. Equation (|3.21| ) can then be written compactly as: 



G G {G-e] 



+ positive powers of G. 



(4.2) 



Dropping the positive powers of G and inverting this equation, we see that there 
are two sheets connected together by a square root cut running between two flnite cut 
points, d and c (see Fig. 4.1). On the physical sheet a further cut (running from b to a), 
corresponding to e^^'^\ connects to further sheets. 




cut of e 



F(h) 



physical sheet 



Fig. 4.1 Sheet and cut structure of G{h) 



To proceed with the solution we relate the function G{h) to the resolvent H{h) by the 
formula: 



h 



et2 - 14 
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Gi{h) G2{h), 



(4.3) 



where Gi{h) and 6*2 (^) are the physical sheet and the only other sheet attached to it by 
the cut of e^^'^^ (see Fig. 4.1). Note that this differs slightly from the form (|3.23| ) derived 
in since there are now an infinite number of inverse powers of G in ( 3.21|) . It can be 
derived by the methods used in 0. 

We now proceed exactly as in 0] . Equation ( ^4.3| ) can be written as 



Fi{h) + F2{h) + H{h) = ln(— (4.4) 

where \iiGi{h) = Fi{h) + H{h). The two sheets Gi{h) and 6*2 (/i) are glued together by 
the square root cut coming from F{h). The combination Fi{h) + F2{h), evaluated on the 
cut of F{h), is twice the constant part of F{h) on the cut (the discontinuous part of F{h) 
is of opposite sign on Fi{h) and F2{h) and is therefore canceled). We thus have the two 
equations 

2^{h) + H{h) =-ln{ ^ ' 



e^2 - t4 (4.5) 
2F{h) + ${h) =-\rvh, 



the first coming from the large limit of the character ( [1.3|) and the second, in view 
of (|3.14[) , being the saddlepoint equation ( p.l2|) . These two equations tell us about the 



behaviour of the function 2F{h) + H{h) on the cuts of F{h) and H{h), respectively. We 
have introduced the notation ^{h) to denote the real part on the cut of F{h), and similarly 
for I/I{h). The principal part integral in ( |3.12| ) is thus denoted in ( [4.5| ) by I/I{h). 

Our object is to now reconstruct the analytic function 2F{h)+H{h) = 2 h\G{h)—H{h) 
from its behaviour on its cuts. To do this we first need to understand the complete structure 
of cuts. We already know the structure of cuts of H{h); it has a logarithmic cut running 
from h = to h = b, corresponding to the portion of the density which is saturated at its 
maximal value of one, and a cut from 6 to a corresponding to the "excited" part of the 
density, where the density is less than one. It thus remains for us to understand the cut 
structure of\nG{h). 

The function G{h) has two cuts on the physical sheet. The first cut, running from b to 
a, corresponds to the cut of e^*^^\ the second cut, running from cut point c to cut point d, 
corresponds to the cut of e^^^"^ (see Fig. 4.1). To see whether \nG{h) has any logarithmic 
cuts, we first notice from (^4.2| ) that G{h) is non zero everywhere in the complex h plane 
except possibly at infinity. Thus for lnG{h) the only finite logarithmic cut points are at 
h = b, defined to be the end of the fiat part of the density (this corresponds to the end of 
the cut of e^^'^^), and possibly the cut point c, defined to be at the end of the cut of e^^'^\ 
The only remaining question is whether this logarithmic cut starting at h = b goes off to 
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infinity or terminates at c. For large h we see from ( [4 .21 ) tfiat \nG{h) = Ine + i-G- 
tfiere is no logaritlimic cut at infinity. We conclude that the cut structure of the function 
In G{h) consists of the cuts corresponding to e^{h) and e^*-'^-* connected together by a 
logarithmic cut, whose cut points are b and c (see Fig. 4.1). 

We now introduce two functions F{h) and H{h) defined by 

F(h) =F(h) -In ^ , and H{h) = H{h) +\n-^. (4.6) 
[h — c) h — h 



They are defined such that they have no logarithmic cuts. The equation (|4.5|) now reads 



2f{h) + H{h) = IniU - et2) + ln(-^^^-^) 



(4.7) 



These two equations define the behaviour of 2F{h) + H{h) on all of its cuts. By standard 
methods we now generate the full analytic function 2F{h) + H{h). There are four cut 
points, a and h defining the cut of H{h), and c and d defining the cut of F{h). Their 
values will be fixed later by boundary conditions. Note that this differs from the fiat space 
limit solved in in that F{h) now has two finite cut points whereas previously it had a 
semi-infinite cut (i.e. d — oo). We generate the full analytic function by performing the 
contour integral 



ds ln(t4-et2)+ln(-(^) 



2F(h)+H(h) =r(h) ^ x . x ■ r x . x 

where for compactness we have defined 



(4.^ 



r{h) = ^/{h - a){h - b){h - c){h - d). (4.9) 

The contour Ch encircles the cut [6, a] , whilst Cp encircles the cut between c and d. They 
are illustrated in Fig. 4.2(a). 



c b — d .j c b +j a 

+i -ve -I 



C 



H 



Fig. 4.2 (a)Contours Cn and Cp. (b)Sign convention for square roots of r{K). 
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The zigzag line between c and b corresponds to the cut of ln( ) , The remaining zig- 
zag hne starting at c corresponds to the remaining logarithmic cut, — ln(/i — c). Expanding 
the contours, catching poles on the way and using the fact that logarithmic cuts have a 
discontinuity of ±z7r, we arrive at 



2F{h) + H{h) =ln^i-^ 
+ r{h) 



ds ^ dsX\n{t^ - et2) ds 



0) 



{h-s)r{s) Ji, {h-s)r{s) J -oo {h - s) r{s) 

(4. 

Fig. 5(b) clarifies the sign convention for the square root r{h) on the real axis above 
and below the cuts. Note that, for the cuts of l/r{h) the signs on the cuts are inverted 
compared to Fig. 5(b), i.e. +i <-> —i. The integrals in ( [4.10| ) are defined to be along the 
upper side. 

To fix the constants a, 6, c and d, we expand ( [4.10| ) for large h and compare the 
resulting power series expansion to that obtained from inverting ( ^.211 ): 

2Fih) + Hih) = 2lne+{^-l)l + 0{^). (4.11) 

The terms of (9 ( -p- ) depend on the as yet unknown positive powers of G in ( [1.2|) . Expanding 
( fl.lOl) for large h and comparing to ([4.11| ) we find the three boundary conditions: 

\n{U~tt2) = -^ (K' + v), (4.12) 



.^ = ^?ii^, (4.13) 

^ - 1 = + ^ L ! + S] , (4.14) 

4 '-4 sn{v, k ) 



where we have defined 



^ =^{a-c){h-d), C = a + b + c + d, 
X Q ( u'\ A I u'\ dn(z;,/c') cn{v,k') 

T=3cn(.,fc)dn(.,fc)--^-^-^ and (4.15) 



2K ^ ' ' 

K and K' are the complete elliptic integrals of the first kind with respective moduli k 



and k' = y/l — k"^. E is the complete elliptic integral of the second kind with modulus /c. 
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E{v, k') is the incomplete elliptic integral of the second kind with argument v and modulus 
k' and sn, cn and dn are the Jacobi Elliptic functions. The nome q is defined by 



-eL 

K 



(4.16) 



Finally k and v are defined in terms of the cut points by 



k 



2 _ {a-h){c-d) 
{a-c){h-d)' 



sn 



a 



a — d 



k'). 



(4.17) 



The final boundary condition will come from the normalisation of the density. 

The density p{h) can be obtained from 2F{h)+H{h) by using the saddle point equation 
2F{h) + ^{h) = — In/i and the fact that the resolvent ( |3.14| ) for the Young tableau can 
be written as H{h) = ^{h) =F inp{h). We thus perform the integrals in ( 4.1U ) and, after 
using the first boundary condition and an identity between elliptic functions^ we obtain 



p{h) = ^ - ^ In 

A IT 



04{^{u - iv),q) 
04{^iu + iv),q) 



(4.18) 



where v is defined in ( |4.17| ) and u is defined by 



u = sn 



-i( i {a-h)(h-d) . 

^^1 {a-h){h-dy >■ 



(4.19) 



Integrating p{h) from 6 to a and equating the answer to 1 — 6 to ensure that the density 
is normalized to 1 (the fiat portion from to 6 gives a contribution 6), leads to the final 
boundary condition: 



-2Kiq 
7r2 e{{Q,q) 



[-E + A-(*'W(„,*')-f^-f2E^)]. 



(4.20) 



The elliptic modulus k and the parameter v are fixed in terms of the coupling constants 
t4 and e, by equations ( [1.12|) and ( fl.20|) . It is natural to use the parameter (3 defined by 
(3^ = so that the elliptic parameters k and v are determined through the two parameters 
^4 and p. Through (f4.17| ) k and v fix the two independent ratios of differences between 
cut points. The scale and position of the cut points are then fixed by equations (|4.13|) and 

(iH) • 



For this and many other relations between Jacobi's elliptic functions and theta functions 



useful for performing the calculations of this section see e.g. [16|,|17]. 
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We now concentrate on the free energy. We choose the definitions for t2, ^4 and e 
introduced in section 2: 

t2 = VXt, U = \, t = , (4.21) 

i 

which imphes that 

t2g = A2(^^j for g>2. (4.22) 

This choice is motivated by several reasons. Firstly we want to have a way of measuring 
the area of the surfaces. The factors of A in the coupling constants have thus been chosen 
in order that the power of A counts the number of squares A making up the surface (see 
eq.( |2.3|) ). Each corner of a vertex corresponds to one of the four corners of a square and 
and so is weighted by A^ . A vertex of 2q legs (where 2g squares meet) is thus weighted with 
A^. The definition of e is such that ^ = as before, so that the elliptic parameters k and 
V are defined in terms of A and {3. Finally, for surfaces of spherical topology, each negative 
curvature defect is compensated for by a precise number of positive curvature defects. In 
particular a negative curvature defect t2q is compensated by g — 2 positive curvature defects 
t2 giving an 

overaU factor of /?2(9-2). The power of t will thus correspond to the number of 
excess positive curvature defects. For a surface of spherical topology which takes precisely 
four extra t2 defects to close the surface we would therefore expect a factor (see (|2.3P ). 

We now use a standard formula from matrix models. Denoting the free energy by 
Z{t, X,(3), Z = e-^'^, we have 

^Z{t,X,l3) = j^{{^TrM')-l). (4.23) 

This identity corresponds to grabbing hold of one of the propagators of a free energy 
diagram. From equation ( |3.15| ) we see that 

{lTrM') = ^ + {h). (4.24) 



where (h) = J dh p{h) h. Using the density ( [1.1 8|) to calculate {h) we obtain, after a long 
calculation, 



dX 4A 



_ k'hn') _ - J-(KS - -)(dn2 + k"cn') 

STTcndn^ ^ 4 27r^ 2^^ ^ 



Stt sn2 27r2sn2 An sn 2n sn tt sn 2 



25) 
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where 

O = -dn^ + 2S2) _ s and $ = snE - T(i^ S - |), (4.26) 

and for compactness we have denoted sn(f , k') by sn and similarly for cn and dn. Every- 
thing inside the square brackets of this expression is a function of the two elliptic parameters 
k and f , in other words a function of A and /3. The ^ outside the square bracket can be 
expressed from ( [4 .131) and (|4.21|) by 

4_0K,ej(JWA (4 27) 

We thus see that the free energy is written as a function of A and (3 times a factor (see 
(p73|)). As discussed above, the factor of corresponds to the four extra t2 defects needed 
to close the surface into the topology of a sphere. 

Since the solution is somewhat complicated, it is essential to check it carefully. As 
a first check we investigate two limiting cases. By setting t2 = e = and t4 = X (or 
equivalently t = P = 1) we correctly recover the solution found in for the model with the 
external matrices A = B = X^J^ where -^Tr[jr'^] equals one for q even and zero otherwise 
(see equations (5.8)-^(5.11) of ||^ ). Taking the other limit and setting e = we correctly 
recover the solution for the flat space limit found in 0. Both these limits are somewhat 
subtle to extract since they correspond to singular points in the boundary conditions, and 
thus have to be derived by a careful taking of limits. We elaborate on this in appendix 
2. where further nontrivial checks are also discussed. 



5. Extraction of physical results 

We are now in a position to derive the physical results announced in section 2. The 
boundary conditions contain sufficient information to find the critical values of the cou- 
plings A and (3 corresponding to the continuum limit (the elliptic parameters k and v do 
not depend on the coupling t which is therefore not fixed by the continuum limit). The 
most direct way of finding the continuum limit, however, is from the expression for the 
density (|4.18|) . The first term of the density, is always positive on the interval (6, a). 



However, for a certain range of values of v, the second term in ( [4.18| ) can be negative. The 



continuum limit corresponds to the point in (g, v) parameter space where the density is no 
longer positive. Specificallyi we need to find the point in [q, v) space where the density 



^ This is very similar to ordinary matrix models. There the point in coupling constant space 
at which the density of eigenvalues becomes flat at its end point corresponds to the continuum 
flmit. 
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becomes exactly flat at the end point a. Taking the derivative w.r.t. h of the density at 
the end point a and setting it equal to gives 



1 



cn(v, k') 



sn(f , k') dn{v, k') 



(5.1) 



where S is defined in (f4.15| ) . This equation defines a line in {q,v) space. Substituting 
( ^T|) into (14.201) , a point on this line can be related to A by 



and /3 by 



7r2 6';(0,g) V cn^v.k') 
1 



1 



A 



(5.2) 



(5.3) 



Equations ( p.l|) , ( W^ ) and ( |5.3| ) define the continuum line in the parameter space (/?, A). 
A plot of this line is shown in Fig. 4.3. 




Fig. 4.3 The critical line in the /3, A plane 



The point /3 = 1 (A = |) corresponds to the model A = B = X^J' mentioned above and 
solved in |jl[]. The point (3 = {X = 1) corresponds to the flat space limit studied in [^. The 
continuum line corresponds to a singularity in the mapping between the coupling constants 
/3 and A and the parameters q and v through which physical quantities are expressed. The 
type of singularity determines the physical regime of the continuum limit. It is easiest 
to investigate the behaviour of the singularity about the limiting points discussed above. 
The limit /3 = 1 was discussed in M , where a standard square root singularity was found, 



18 



corresponding to the pure gravity regime. We thus concentrate on the hmit /3 ^ and 
A ^ 1. 

This hmit corresponds to q 1. We can study it most simply by expressing ah the 
dependence on the nome q through the dual parameter q' , defined through 

= e""""^ = e~i^. (5.4) 

We then discard all exponentially suppressed terms, i.e. all powers of q' . This is not the 

2 

same as setting q' = since we keep terms of order Ing', in particular K = — (1 + 
0{q')). For q close to one this is a good approximationH. In this limit the boundary 
conditions fixing v and the nome q are given by 



2zv 

-e ^ 



(1 H ) cos2t; + - Iz H l) sm2v 

TV IT TT 



(5.5) 



where we have defined z — — ln(A(l — /3^)). Setting the derivative (w.r.t. v) of the first 
of the above equations to zero gives the critical line about the limiting point where the 
surface fiattens: 

z 

tanfc = — . (5-6) 

TT 

This could also have been obtained by discarding all powers of q' in . To investigate 
the type of singularity on the critical line we expand the r.h.s. of the first equation in ( |5.5| ) 
about the critical value of v: 

^ ■■ao + a2{v-v,f + 0{iv-v,f). (5.7) 



1-/32 

The coefficient gq is given by evaluating (|5.5| ) at Vc- The coefficient 02 is the second deriva- 
tive (w.r.t v) of the first equation of (|5.5|), also evaluated at Vc- It is easy to demonstrate 
that a2 cannot be zero for positive, real coupling coefficients. We thus find that the singu- 
larity on the critical line is a square root corresponding to pure, two dimensional gravity. 
A more rigorous study would analyse the singularity of boundary conditions ( [4.12| ) and 
( [1.20|) at a generic point along the critical line. We have examined this numerically and 
find a square root singularity for the full length of the critical line. We can thus conclude 
that for all finite (3 we are in the regime of pure, two dimensional gravity. 



^ Numerically the computation of elliptic functions is performed by writing them as a power 
series in either q or q' , whichever is smaller. The power series converge so rapidly that, even for 
a generic point, the first few terms are sufficient to give a high degree of accuracy. 
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We now proceed to extract the physical mechanism of flattening which requires the 
hmit /3 ^ (corresponding to f ^ 0). Expanding out ( |5.5| ) to lowest order in f , z and fi'^ 
we flnd 

V — (x — a/x^ — l) and q — X. (5-8) 
where we have introduced a scaling parameter x defined by 

x = l + — ^ with /x = Ae-A and = 1 - ^ + 0(/3^). (5.9) 

Note that this is a very natural, dimensionless scaling parameter; /i is the continuum cos- 
mological constant with dimension of inverse area and (3 controls the number of curvature 
defects per unit area (and thus also has dimension of inverse area). 

Expressing the elliptic functions in ( [4.25|) in terms of the dual nome q' , dropping all 
powers of q' and keeping only the leading order contribution in powers of v, we arrive at 
the free energy about the flattening transition point: 

Z{t, A, /3) = 3^ [x' - \x' + y - ^ - - If"] ■ (5.10) 

This is the continuum scaling function corresponding to the flattening transition. For x 
close to 1, in other words for cosmo logical constant much less than /3, there is the standard 

5 

^2 singularity characteristic of pure gravity. For x very large, corresponding to very small 
13 the free energy has the singularity //~^. From (|5.10| ) we can also Laplace-transform to 
the free energy for fixed area 



Z(t,A,/3) 



/48 8\ , , /24 , . 



t4^-fe^^'^A-i, forA^oo 
t^A, for A ^ 



TT 

with y = -j=PA, 

^2 (5.11) 



where Iq and Ii are modified Bessel functions. The limiting behaviours for large and small 
areas are also shown. For surfaces whose area. A, is large compared to (3 the exponent 
— I signals that we are in the regime of pure 2D gravity. For surfaces of very small area 
A there is a linear dependence on the area. The function Z{A,P,t) smoothly interpolates 
between these two limiting behaviours. There is no phase transition separating flat space 
from pure gravity. 
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6. Conclusions and open issues 

In conclusion, our result strongly suggests that gravity (|1.1|) does not seem to exist 
as a continuum theory. The discretisation of the term by the same principle that 

consistently defines the path integral ( |l.lj ) for (3o = oo (i.e. ordinary 2D gravity) does not 
lead to any new fixpoints, indicating the non-perturbative irrelevance of the dimensionful 



higher derivative terms of ( |1 . 1|) . Of course one might object that our particular (but 
generic!) discretisation is too naive. Indeed, it would be gratifying to prove the universality 
of our result by choosing different weights ( p.2|) for the curvature terms in our model. In 
principle, as we hope we have convinced the reader, this can be checked with the help of 
our formalism. It is not entirely excluded that fine-tuning the weights (|2.2|) in a subtle 
way might result in a smooth phase of gravity (for some (non-rigorous, as we believe) 
arguments in favor of this hypothesis see |]T^). However, it is equally evident that any 
further calculations will be very tedious unless it proves possible to simplify our method 
in a significant fashion. We should also point out that the fact that the curvature in 
the present type of discrete models is always necessarily "quantised" in units of tt could 
be an a priori stumbling block for reaching a smooth phase of gravity. Could one find 
discrete models of gravity that allow for soft, slowly varying curvature fiuctuations over 
many lattice spacings? 

The present work presents the first rigorous step towards a full understanding of these 
issues. It would be very interesting to gain deeper insights by simplifying and extending 
our novel discrete approach. 

One might also attempt to describe a model like ours directly in the continuum. One 
could start with the conformal metric of a fiat surface with localised curvature defects. It 
can be represented locally as gab = Sab e'^^^^ with 

M 

^{z) = Y, Rj H^-^jf. (6.1) 

where Rj = —1,1. Symbolically, the partition function might be written as 

M ^ 

Here we introduced the fugacity of curvature defects f3 instead of the explicit i2^-term in 
the action. It serves the same purpose: for — we arrive at the completely fiat metric, 
whereas for /3 ~ 1 the system should show the behaviour of pure quantum gravity, at least 
in the infrared domain. We retained the notation /3 to denote the parameter playing a role 
similar to the R^ coupling in the above discrete model. 
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This formulation resembles a little bit the two-dimensional Coulomb gas problem. 
However, the measure of integration d[zi, ...,zm] of the positions of the curvature defects 
is a complicated object: it should take into account the topology of the surface and the 
existence of zero modes (the action does not depend on some directions in the space of the 
Zi). It would be very interesting to make this direct continuum formulation more precise. 

Another interesting issue is the role of exponential corrections appearing due to the 
structure of elliptic functions. In fact all physical quantities, such as Z{t, A, (3) (see (|2.3|) )), 
contain exponentially small terms in the limit A ^ Ac and /3 0, thus leading to an 
essential singularity at /9 = 0, A = 1. These terms are precisely the powers of q' that we 
discarded to arrive at the scaling function ( p.lUD of section 5. One can obtain the first 
correction of this type in e.g. the free energy /(/3) per unit area in the ther mo dynamical 
limit A = Ac: 

Z(t,A/3) ~ ^e^"(^)^ A-i with f((3) = ^(3\{l + ...)+e-^{4 + ...)] (6.3) 
p2 V2 L J 

where f{(3) = limA^oo ^ ln2^(t, A, /3) = Ac and the dots denote terms of order and 
higher. 

These exponential terms are likely to be lattice artifacts. They emerge even in the 
simplest calculation for the flat closed quadrangulation with four positive curvature defects, 
where they appear as discrete corrections to the approximation of elliptic sums by integrals 
close to the continuum limit. 

On the other hand, formula ( |S.3| ) corresponds to the critical free energy as a function 
of the curvature fugacity f3 (i.e. we have already taken the continuum limit). It is possible 
that the exponential terms might be corrections relevant for the statistical mechanics of 
random lattices at long distances (of order ■^) rather than for continuous 2D-gravity. 

However, the most interesting and rewarding extension of our methods and ideas would 
be their generalisation to allow for the coupling of matter to discrete two-dimensional 
gravity. It would lead to an entirely novel approach to two-dimensional physics, unifying 
theories in the absence and presence of two-dimensional quantum gravity. 
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Appendix 1. General case: Large limit of Schur characters 

The two equations, (|3.25| ) , in which there are three non-zero couplings, ^4 and tQ, 
can be reduced to a single integral equation for the density. The derivation we present below 
can be generalised in an obvious manner to an arbitrary finite set of non-zero couplings. 
Other models with different saddlepoint equations can also be treated with only minor 
modifications. 

The initial steps are very similar to the start of section 4, to which we refer the reader 
for more details. We start from the equation 

h — 1 = + + + positive powers of G. (1.1) 

and deduce (by identical reasoning to section 4) that the physical sheet of the function G{h) 
and all the sheets attached to it by the cut of e^^'^^ have a cubic structure. We further 
deduce that the function In G{h) has a logarithmic singularity starting on the physical 
sheet at the point h = b, and, by looking at the large h limit of ( p. . 1| ) (inverting ( |1 . 1| ) to 
lowest order in gives \iiG{h) = ^ In(^) + 0{h~^) ) deduce that the logarithmic cut goes 
off to infinity. We now study the combination 2lnG{h) — H{h) where H{h) is defined in 
( [1.6|) . We find it has the analytic structure shown in Fig. Al 




Fig. Al 2\n.G{h) — H[h) in the complex h plane 



where the zig-zag line, starting at /i = 6 on the physical sheet, represents the logarithmic 
cut. It also has a saddlepoint equation (trivially derived from ( |3.25| ) and ( ^^ ) ) of 

lk[2\nG{h)-H{h)] = -\n{h-h). (1.2) 
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We now unwind the three sheets of Fig. Al by parametrising /i as a cubic polynomial of a 
new parameter z 

h = z^-az + '-i. (1.3) 

The coefficients a and 7 are such that the cut points of /i as a function of z are identical 
to the cubic cut points of Fig. Al. The three sheets of Fig. Al (defined in the complex h 
plane) then unwind to become a single sheet with three cuts (plus a logarithmic cut) in 
the complex z plane (see Fig. A2). 




Fig. A2 2\nG{h{z)) — H{h{z)) in the complex z plane 



Defining 



J^{z) = 2\nG{h{z)) - H{h{z)) - 21n(- 



'■6 



hi 



:i.4) 



we see that J^{z) is an analytic function in the complex z plane with three cuts corre- 
sponding to the images of the cut of the density. Furthermore, we see from the large h 
limits of \rvG{h) and H{h), given by \nG{h) = i In(^) + 0{h-^) and H{h) = C(i), that 
the function J-'{z) behaves as 0{^) for large z. The function J-'{z) can thus be generated 
entirely from the discontinuities of its cuts: 



61 z- s 



-,,P(M£)) 

62 z-s 



(1.5) 



We now use the identity {zi — z'i){zi — 2^2) (^1 ~ ^3) = h — h' = {zi — z'i){z2 — z'i){z3 — z'l), 
where the Zi are the three roots of (|1.3 ) and the z[ are the three roots of (|1.3| ) with h set 
to h' . Along with equation ( |1.4| ) , this leads to a saddle point equation in the complex z 
plane: 



'bi 



ds p{h{s)) 



1 



1 



z — s 



Z2[Z 



zsiz) 



ln[ 



{z-h] 



{z-b2){z-bs) 



-Inte, (1.6) 
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where z, and 2:3 are the three roots of ( |1.3|) . z^ and zj, are thus functions of the first 



2 



root z. With the parametrisation of (|1.3| ) they would be given by 2:2(3) (-2) = f i y — 

The generahsation to more than three weights is obvious. However, the functions p 

that satisfy such an equation, even for the case of just three weights, involve functions 

more general than elliptic and probably more general even than hyperelliptic. 

As a final comment, notice that by unwinding the analytic structure of the function 

hiG{K) using the mapping (|1.3|) as in the discussion above we arrive at a closed expression 

for the function \^G{K) = F{h) + H{h): 

lnGW= m A. (1.7) 



where z is related to h by (|1.3| ). Given a representation, specified in the large limit 
by the density p{h) (which then defines H{h) through ( |3.14|) ), we thus have a closed 
expression for the logarithmic derivative of the character, F{h), ( |3.13| ). The coupling Iq, 
is incorporated directly into the expression (|1.7|), the two remaining couplings, t2 and t4. 



determine, through (|3.21| ), the coefficients a and 7 of the mapping (O). It is again trivial 
to generalize this result to more than three weights. 



Appendix 2. Verification of the solution 

2.1. Limit (3-^1 

This corresponds to the model with A = B = X^J' studied in . Setting t2 = e = 
and t4 = X (or equivalently t = f3 = 1) the first boundary condition, (|4.12| ), implies 
immediately that K'^oo^k = 0^c = d. The values of E, K, E{v, k'), sn{v, k'), S 
and T are then given by 



E = K=^, E{v,k')=sn{v,k') = l, E = 2 and T = -./^-J^,. (2.1) 
2 \a — o\a — c 



where we have used the definition of v in terms of the cut points, ( [1.17|) , along with 



standard identities cn^ = 1 — sn^ and dn^ = 1 — /c'^sn^ to find T. After using ( [4.13|) to 
eliminate the theta functions in ( [4.20|) , and also using ( [1.14|) , we find that c = d = 0. In 



other words the cut of F(h) disappears. We also find the first condition fixing a and b: 

a -3^-1 = where a=( ^+^ )' and ^=( ^^^ )'. (2.2) 
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The hmit of the third boundary condition, ([4.13D , is more subtle. We shift the argument 
of the theta function using the identity: 6i{z, q) = qe~^^^6i{^^^^ — z, q) and carefully take 
the limit -^{K' — f ) — » sinh~^ \f^^ ^° final boundary condition 

SAV^- (7 + 1 = 0. (2.3) 

These were the conditions found in |jl[] leading to Ac = |. The density ( ^l.ltjj ) can likewise 
be reduced down to that found in [Q. 

2.2. Limit /3 ^ i.e. e ^ 

This limit recovers the flatspace solution found in 0. In this limit d — oo and the 
density ( |4.18[ ) trivially reduces to that found in 0. To see that the boundary conditions 
give the correct limit and also to check perturbatively that the next few powers in /3 
correctly correspond to negative curvature insertions we expanded the results in powers 
of V. Expansions of elliptic functions can be found in most mathematical tables. A 
comprehensive list is found in |TH]. Below we give the expansion formula for the theta 
function which is less common 

M|-).-„_-(:,.._3|K,o,A ,2.4, 

Expanding out (|4.12|) and ( [4.2q ) leads to 



v = ^P^ + 0{P^) and q = U-^-^+0{P^). (2.5) 

For /3 = we recover the first boundary condition of . Using the identities (a — c) = 
^"cn" ^ ^ ~ along with ( p.5|) and (|4.13[) and working to lowest order we find the 
second boundary condition of 0: 

^ Va^. (2.6) 



A careful expansion of ( [4.14] ) then leads to the final boundary condition found in 



+2 

a^l + ^{K^ - EK). (2.7) 
A final check of the solution can be performed by expanding out the expression for the 



free energy ( [1.25|) in powers of v and hence in powers of /3. To lowest order we correctly 
recovered the free energy found in 0. Expanding out to the first non-zero order in P we 
retrieved the order contribution. The result matched exactly with the calculation of a 
single negative curvature insertion of t^, calculated in M. 
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